Factorization and entanglement in general XYZ spin arrays in non-uniform transverse 
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We determine the conditions for the existence of a pair of degenerate parity breaking separable 
eigenstates in general arrays of arbitrary spins connected through XYZ couplings of arbitrary range 
and placed in a transverse field, not necessarily uniform. Sufficient conditions under which they are 
ground states are also provided, ft is then shown that in finite chains, the associated definite parity 
states, which represent the actual ground state in the immediate vicinity of separability, can exhibit 
entanglement between any two spins regardless of the coupling range or separation, with the reduced 
state of any two subsystems equivalent to that of pair of qubits in an entangled mixed state. The 
corresponding concurrences and negativities are exactly determined. The same properties persist in 
the mixture of both definite parity states. These effects become specially relevant in systems close 
to the XXZ limit. The possibility of field induced alternating separable solutions with controllable 
entanglement side limits is also discussed. Illustrative numerical results for the negativity between 
the first and the j th spin in an open spin s chain for different values of s and j are as well provided. 

PACS numbers: 03.67.Mn, 03.65.Ud, 75.10.Jm 
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Quantum entanglement constitutes one of the most 
fundamental, complex and counter-intuitive aspects of 
quantum mechanics. It is an essential resource in quan- 
tum information theory [lj, playing a key role in quan- 
tum teleportation Q and computation [H, d, 0| • It also 
provides a deeper understanding of quantum correlations 
in many-body systems In particular, a great effort 
has been devoted in recent years to analyze entanglement 
and its connection with critical phenomena in spin chains 
[1, H, 0, HI- Studies of finite chains, of most interest for 
quantum information applications, are presently also mo- 
tivated by the possibility of their controllable simulation 
through quantum devices [H, [l(| • 

A remarkable feature of interacting spin chains is the 
possibility of exhibiting exactly separable ground states 
(GS) for special values of the external magnetic field, first 
discovered in [Il|,[l2j m a 1^ XYZ chain with first neigh- 
bor coupling. It was recently inve stig ated in more general 
arrays under uniform fields [H, [MEl, El, E3, [ll, Ell, 
with a completely general method for determining sepa- 
rability introduced in (Tgj |. Another remarkable related 
aspect is the fact that in the immediate vicinity of these 
separability points (SP) the entanglement between two 
spins can reach infinite range [H, [TtJ • In [I?) we have 
shown that the SP in finite cyclic spin 1/2 arrays in a 
transverse field corresponds actually to a GS transition 
between opposite parity states (the last level crossing for 
increasing field) , with the entanglement between any two 
spins reaching there finite side limits irrespective of the 
coupling range. In a small chain, this SP plays then the 
role of a "quantum critical point". In contrast, the en- 
tanglement range remains typically finite and low at the 
conventional phase transition Q. 

The aim of this work is to generalize previous results 
to XYZ arrays of arbitrary spins and geometry in a gen- 
eral transverse field, not necessarily uniform. Moreover, 
we will also determine the exact side limits of the entan- 



glement between any two subsystems (including those for 
the block entropy and those for any two spins or group of 
spins) at the SP analytically, for any spin value. A non- 
uniform field will be shown to allow exact separability 
with infinite entanglement range in its vicinity in quite 
diverse systems (such as open or non- uniform chains), 
including the possibility of field induced alternating sep- 
arable solutions along separability curves, with control- 
lable entanglement side limits. Illustrative results for the 
negativity between the first and j th spin in an open spin 
s chain as a function of field and separation are as well 
presented, for different spin values. 

We consider n spins Sj (which can be regarded as qu- 
dits of dimension d% = 2s, + 1 > 2) not necessarily equal, 
interacting through XYZ couplings of arbitrary range in 
the presence of a transverse external field b 1 , not neces- 
sarily uniform. The Hamiltonian reads 



H 



iQ^ + ^+t,^), (1) 



and commutes with the global S z parity or phase-flip 
P z = exp[i7r^" =1 (sf + Si)] for any values of b l , or 
Si. Self-energy terms (i — j), non-trivial for S, > 1, 
are for instance present in recent coupled cavity based 
simulations of arbitrary spin XXZ models [l(| and will 
be allowed if Sj > 1. 

We now seek the conditions for which such system will 
possess a separable parity breaking eigenstate of the form 



I©} 



tf =1 exp[i0 iS V]|O 

fe=0 



^fsin*!^)], 



(2) 
(3) 



where sf|fej) = (k — Si)\ki) and e l9iSi |0j) is a rotated 
minimum spin state (coherent state [H|). The choice of 
y as rotation axis does not pose a loss of generality as 
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any state e l ^ i ' &i \Qi) corresponds to a suitable complex 
Bi in © [21]. Replacing sf in ||T]) by e" 



Definite parity eigenstates of H in the subspace gener- 
ated by the states | ± O) can be constructed as 



z,x 

i.e., V 



s\ ,x cos Bi ± sf' z sin^i, — > sj , the equation 



ff|0) = £ e |©>, i.e.^elO) = £ e |0> with |0> = ®? =1 |0i) 
and i?e = e~ l ^ SiS i }{ e l Y.i s i s i ; leads to the equations 



iC sinfi sm( 



6 l sin 0, 



(sj — i5y ) (v l J cos 0^ sin Bj — v l J sin (9, cos 



which determine, for instance, the values of Vy J and V in 
terms of v % J , v l J , and 9i . The energy is then given by 



E e = 



E 



s t [b % cosOi + | Y^(sj - |<%)(^ sin^ sinfy 



cos fl ( cos + +<+<)]• 



(6) 



For a II? spin s cyclic chain with first neighbor couplings 
( w u = u M^«>i±i) m a uniform field (6 1 = b) we recover the 
original GS separability conditions of ref. [II] for both 
the ferromagnetic (v^ > 0, Bi — B) and antiferromagnetic 
(Vp < 0, Bi = cases. Eqs. fll}-© are however 

completely general and actually hold also for complex 
values of Bi, v l J and b l : If satisfied Vi,j, will have a 
separable eigenstate (J2J) with eigenvalue If sin0j ^ 
for some «, this eigenvalue is degenerate: |0) will break 
parity symmetry and therefore, the partner state 



e)=P,|0) = ®g =1 ejtp[-*e i a?]|O < 



(7) 



will be an exact eigenstate of if as well, with the same en- 
ergy ([6]) . The points in parameter space where the states 
| ± 6) become exact eigenstates correspond necessarily 
to the crossing of at least two opposite parity levels. 

For real Bi, Eq. is just the stationary condition for 
the energy ([6]) at fixed b 1 , t$. The state ([2]) can thus 
be regarded as a mean field trial state, with Eq. ([5]) the 
associated self-consistent equation. Eq. (U), which is spin 
independent (at fixed v% ) , ensures that it becomes an 
exact eigenstate by canceling the residual one and two- 
site matrix elements connecting |0) with the remaining 
states. Moreover, if Bi S (0, ir)\/i and 



|t#'l<»£ Vi,j : 



(8) 



we can ensure that jtO) will be ground states of H : In the 
standard basis formed by the states \h}}, the terms 

in H depending on {sf} are diagonal whereas the rest 
lead to real non-positive off-diagonal matrix elements, as 



E^/iEfsfSj = T,v=±< j {4 s j " + s * s j)> where 4 



M "i "J 

s* ± is) and v% = \(v% ± t$) > by Eq 



. Hence, 

(H) can be minimized by a state with all coefficients real 
and of the same sign in this basis (different signs will not 
decrease (H)), which then, cannot be orthogonal to |0) 
(Eq. |3])). With suitable phases for Bi, | ± 0) can also be 
GS in other cases: A n rotation around the z axis at site 
i leads to Bi — > —Bi and v % J — > — for z 7^ j. 



|0±) = 



|6)±|-9) 
v/2(l±Oe) : 



(4) 



O e = (-0|0)=nr=icos : 



2si 



(9) 
(10) 



which satisfy P 2 |0 ± ) = ±|0 ± ), (©""I©"') = S vv> . Here 
we have set Bi real Vi, since by local rotations around 
the z axis we can always choose yi in the direction of (pi 
(and hence Bi real) in the final state |0). Moreover, we 
may also set \9i\ < tt/2 (and hence Oe > 0) since a local 
rotation of 7r around the x axis leads to 9i — ► 7r — The 
overlap (|10p will play an important role in the following. 

When the degeneracy at the SP is indeed 2, the states 
© (rather than @) are i/ie actual side limits at the SP 
of the corresponding non-degenerate (and hence definite 
parity) exact eigenstates of H. For small variations 5b l , 
the degeneracy will be broken if Oq 7^ 0, with an energy 
gap given by AE « 5b l AMi, where 



am, = <e-|«? !©-)-<© 



l©+) = 



2s, ; sin 2 8 t O e 
cos^l-O 2 ,) ' 



(In contrast, (±0|sf|±0) = -s.cos^). When (©*) are 
GS, a GS parity transition |0~) — > |@ + ), characterized 
by a magnetization step AM — AMi, will then take 
place at the SP if all or some of the fields are increased 
across the factorizing values ([5]). If AE or AM can be 
resolved or measured, the realization of the states ([9]) 
is then ensured. Their magnitude is governed by the 
overlap IT0|) . appreciable in small systems (if Bi ^ tt/2) as 
well as in finite systems with small angles 8f « Si /n, such 
that O e ~ e -^ 8iSi / n . This implies (Eq. gj) systems 
close to the HZ limit = uj/). In this limit {B l -> 0), 
AM 1, with |©+) |0) and |0-) cx EiV^i^l 1 *) 
(weighted VF-type state), where |lj) = <8"=i|(^i)j')- 

Entanglement of definite parity states. In contrast with 
|± 0), the states ([9]) are entangled. If sin#j 7^ V i the 
Schmidt number for any global bipartition (A, A) is 2 and 
the Schmidt decomposition is 



|e ± ) = 
v% = 



(1+vOaK1±vO n) 
2(l±Oe) 



(11) 

Oa = (-©a|©a> , (12) 



where |0^), |©^) denote the analogous normalized defi- 
nite parity states for the subsystems A, A, with v = ±, 
= Oe and + = 1. Hence, I©"* 1 ) can be 
effectively considered as two qubit states with respect to 



any bipartition (A, A), with 1©^), |@^) representing the 
orthogonal states of each qubit. Accordingly, the reduced 
density matrix of subsystem A in the state {Q^) is 
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P A Z + \ Q A)(®A\+Pi-\®1)(®A\ 



(13) 



The entanglement between A and its complement A 
can be measured through the global concurrence (square 
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root of the tangle [22$ ) C AA — \J1{\ — trp 2 ^), which 
for a rank 2 density is just an increasing function of 
the entanglement entropy E AA = — tr pa log 2 pa, with 
= E AA = (1) for a separable (Bell) state. In the 
states (fTTj) we then obtain 



support will lie in the subspace spanned by the four states 
{\®b)\ Q c) > v > u ' = ±}> such that 



{l-0 A ){l-0\) 



Pa+1bc+ 




l±O f 



(14) 







+ ^BC 





P^lBC+ P^- a J3C 

P A - a BC Pa-1~bc- 





\ 



p% a ic 




Pa+1bc- J 



where " 



These values represent the side limits of C AA at the SP. 
For O e > 0, C AA > C+ A , with C AA = 1 if A = O a . 
Note that |0 ± ) are simultaneous Bell states for (A, A) 
only if O a = C?4 = (GHZ limit of |6 ± )). Increasing 
overlaps will in general decrease the global entanglement. 

At the SP, the entanglement entropy of a block of L 
spins in a ID first neighbor spin 1/2 XY chain in a con- 
stant field was found in [23| to be Sl = — trpL In pl = hi 2 
(i.e., C L i = E L i = 1) in the thermodynamic limit, in 
agreement with Eq. (|14[) for vanishing overlaps. Eq. 
extends this result to general finite chains, lead- 
to a slightly smaller value: For small Oa, O a , 



1bc» 



and 



(l+uO B )(l±vOc) 
2(1±ObOc) 



l BC 



Qbc+Qbc- 



BC+ 



1. p A will be entangled if its 
partial transpose has a negative eigenvalue [281 ] . a con- 
dition here equivalent to a positive mixed state con- 



currence 



C 



BC 



Max[C|,Cf ,0], where C, 



± _ 



2\p A \a BC — p A _„a B v c ] represent parallel [y = +) or an- 
tiparallel {y = —) concurrences, i.e. driven by \@ A ) or 
|8^) in Eq. (TJ). This leads to Eq. JT5J, with C+ c (C£ c ) 



parallel (antiparallel). The ensuing negativity, given here 
by minus the negative eigenvalue of the partial transpose 
|,iS is then given by Eq. (fT6|) . 
For B = A, C = A (Og^ -> 1), Eq. (TJJ| reduces 



(Pa) 



C± « 1 - i(CU ± O^) 2 and Sf « In 2 - i(O z ± £ ) 2 to d> with ^ = For a P air of s P ins ^ J> 



AA 

(with L = (^)* in the s = 1/2 XY chain). 

Pair-wise and subsystem entanglement. On the other 
hand, the entanglement of a subsystem is enabled by non- 
zero overlaps. A remarkable feature of the states © is 
that any two spins or disjoint subsystems B, C will also 
be entangled if the complementary overlap O b+c is non- 
zero and O b < 1, O c < 1. Moreover, this entanglement 
can be characterized by the concurrence 



Ob = cos 1 9i, Oc — cos ls i 9j and the result of [171 is 
recovered from (fT5|) if s; = i and 6*,; = 6* Vi. We finally 
note that if B = O c , N+ c = C^ c /2, as in the case 
of a global partition. In general, however, there is no 
proportionality between N BC and C BC . 

The concurrences (|15[) fulfill the monogamy inequali- 
ties [31( C B C+D > C BC + C BD for any three disjoint 
subsystems B,C,D. We actually obtain here 



^BC — 



^{i-a B ){i-ol)o 



B+C 



l±O e 

or equivalently, the negativity [H, [25[ , 



(15) 



Cbc 



c 2 — c 

^BD — ° 



B,C+D 



(l-0 2 r: )(l-0 2 n ) 



i-o%o 



7- 



(17) 



Let us also remark that subsystem entanglement per- 
sists, though attenuated, in the uniform mixture 



N 



BC 



^(P%) 2 + i^c) 2 /(>B+U-pU^ (16) 



P° = i(|e+>(e+| + |e->(e-|), 



(18) 



where A = B + C. While the concurrence of an ar- 
bitrary mixed state pa (which can be defined through 
the convex roof extension of the pure state definition 
[26l ]) is not directly computable in general (the excep- 
tion being the case of two qubits [27]]), the negativity 
Nbc = ^ [Tr|p^f I — 1], where p A denotes partial trans- 
pose [HI, can always be calculated [2l|, being just the 
absolute value of the sum of the negative eigenvalues of 
p A . Eq. (fT6"]) represents then the side-limits of Nbc a t 
the SP. 

Proof: For A = B + C, we first note that if O a = 0, 



Eq. (QD becomes p\ = ±(\Q A )(Q A 



&A)(-e A \ 



i.e., p A coincident and separable (convex combination of 
product densities [13]). Entanglement between B and 
C can then only arise if O b+c 7^ 0. Next, using similar 
Schmidt decompositions (fTTj) of the states |0^), Eq. JT5|) 
can also be considered as an effective two-qubit mixed 
state with respect to any bipartition (B, C) of A: Its 



which differs from |(|0)<0| + |— 0><— 0|) if O e ^ 
and represents the T — » + limit of the thermal state 
p cx e~ m at the SP when | ± 6) are GS (and the GS 
degeneracy there is 2). Replacing p AlJ by h(p A " +Pa") m 
()13|) . we find now antiparallel global and subsystem con- 
currences, given for any disjoint subsystems B, C by 

C°bc = h( C Bc - C%c) = C BC Oe/(l + O e ) , (19) 

i.e., half the parity splitting of Cbc- Eq. (fT9"j) remains 
valid for a global bipartition (B — A, C — A). The 
ensuing negativity can be similarly calculated. 

The order of magnitude of subsystem concurrences is 
governed by the complementary overlap O B+c - For small 
subsystems (like a pair of spins) in a large system, C BC 
will be appreciable just for sufficiently small angles in 
the complementary system, i.e., 9 2 w Si/n, such that 



a 



B+C 



- J2 ie A s i S i/ r - 



remains finite. This leads again 



to systems with small XY anisotropy. 
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Uniform Solution. Let us now examine the possibility 
of a common angle 9i — 9 V i. Eq. © leads then to 



COS 3 



(20) 



implying a fixed ratio x = ( v y ~ v l J )/( v x ~ V V) — c °s 2 
for all pairs with v l J ^ v l J , and an isotropic coupling 
v y = v x ^ = v% i ■ ^ subset of isotropic couplings 
will not spoil this eigenstate [32|. Eq. (O implies then V 
arbitrary if 9 — or it (XXZ case v l J = v l J) or otherwise 



The energy © becomes 



' ) ( S j 2 ) 



(21) 



(22) 



A general field allows then a uniform separable eigenstate 
(a global coherent state) in cyclic as well as open chains 
with arbitrary spins Sj in any dimension if (|20p holds 
For instance, in an open II? spin s chain with 
first neighbor couplings v l J = v^Ajii, Eq. (|2"Tj) yields 
b l = b s = 2syJ{v y — v z ) (v x — v z ) at inner sites but b 1 — 
b n = \b s at the borders. 

Eqs. (j20 |) -([22 |) are actually valid for general complex 
9, but real fields imply cos 9 real (x > 0). The case 
cos 2 9 > 1 (imaginary 9) corresponds to a rotation around 
the x axis but can be recast as a rotation around the y 
axis by a global rotation around the z axis. Hence, we 
may set cos 2 9 G [0, 1]. | ± 6) will then be GS when Eq. 
© holds. 

The concurrence (fT5)) becomes, setting cos 2 9 = x, 



^bc — 



^(1 _ x 2 S B )(l _ x 2S c } x S-(Sb+Sc) 

i±x s 



(23) 



where Sb — J2ieB s i is ^ ne subsystem total spin and 
5 = s i t ne total spin. It is independent of separation 
and coupling range, depending solely on x and the ra- 
tios S B /S, S c /S. If x = 1 - 5/2S, with 6 > and finite, 
X s ~ e~ s / 2 remains finite for large 5. Eq. (j2"3"]l leads then 
to 0(1/ VS) and 0(1/ S) global and subsystems concur- 
rences for small Sa, Sb and Sc- 



C 

c 



± _ 

AA 
± 



5a<5 Vl - e~ s 
S l±e- 5 /2 

g VSB~S^e- s / 2 
S 1 ± e" 5 / 2 ' 



(24) 
(25) 



On the other hand, for Sa — |5, C AA - — 1 whereas 



AA 



tanh j5. Thus, while for large 5 both C AA rapidly 
approach 1 as Sa increases, for small S (XXZ limit) this 
occurs just for C AA and Sa close to 5/2 (here |6 + ) — > |0) 
but |6~) approaches the VF-type state cx y^sjlj)). 

Alternating solution and controllable entanglement at 
the SP. Among other possibilities allowed by Eqs. ((I])-©, 



let us examine that of a field induced two-angle solution 
in a ID chain (cyclic or open) of spin s with first neighbor 
XY couplings (v l J = <5jj±iw M , with v z = 0). We assume 
X = v y /v x € [0,1]. A separable eigenstate with #2i = 9 ei 
= 9 a is feasible if there is an alternating field b 2t = 
1 = b Q in inner sites satisfying (Eqs. ©-©) 



b P , b 



2(- 



bJ)n 



(2s) 2 v x v y 



(26) 



This leads to a transverse separability curve. The ensuing 
angles satisfy cosfloCOS^e =v y /v x and are given by 



cos 



bl + (2sv x y 



(27) 



being field dependent. For b e = b Q we recover the previ- 
ous uniform solution (b s — 2s^Jv x v y ). In an open chain 
we should just add, according to Eq. ©, the border cor- 
rections b 1 = 2& ) b n = \b an . The states | ± 0) will then 
be GS setting 6* , e > when v x > and 9 a > 0, 9 e < 
in the antiferromagnetic case < (for even n if chain 
is cyclic, to avoid frustration). 

The definite parity states |0 ± ) will again lead to in- 
finite entanglement range, but with three different field 
dependent (and hence controllable) pairwise concurrences 
between any two spins (Eq. (fl~5| for B = i, C = j): even- 



even, odd-odd and even-odd, satisfying Cf e = y C^bC^k, 
with Cf Q > C± e > C± if \b a \ < \b e \. Hence, Cf Q can 
be made larger than C^ e despite the absence of odd-odd 
direct coupling. For sufficiently large b e , cos9 e w 1 but 
cos 9 w x : just odd-odd pairs will be appreciably entan- 
gled in this limit at the SP. 

Application. As illustration, we first depict in fig. [1] 
full exact results for the GS negativities N\j between 
spins 1 and j in a small open chain of uniform spin s 
with first neighbor XY couplings in a uniform transverse 
field b % — b for i — 2, . . . , n — 1, with the border correc- 
tions b 1 = b n = ifr. For x — v y/ v x £ (Q> 1) this chain 
will then exhibit an exact factorizing field b s = 2sv xy /x 
where separable parity breaking states with uniform an- 
gle cos 9 — y/x will become exact GS if v x > (if 
v x < 0, 9i = (—l) l 9 instead in the GS). We have set 
X = 1 — 5/ (2ns), such that the side limits of the negativ- 
ity at b s are roughly independent of s and n. It is first 
seen that the ensuing behavior of the N\j in terms of the 
scaled field b/b s is quite similar for the three spin val- 
ues considered (s — 1/2, 1 and 3/2, the latter involving 
a diagonalization in a basis of 65536 states for n = 8). 
The GS exhibits ns parity transitions as the field is in- 
creased from + to b s , with the last transition at b s . As 
the latter is approached, it is verified that the pairwise 
entanglement range increases, with all negativities ap- 
proaching the common side limits (|16jl . distinct at each 

side, given here by N£ = \Ct ~- 



and NZ 



*3 

5 2 e- 5 / 2 



2n(l+e- 5 / 2 ) 



(Eq. ((251)) 



An interval of full 



•n ~ 4p"_ 4n 2 (l-e- 5 /2)2 ■ 

range pairwise entanglement around b s is then originated, 
which involves on the left side essentially the last state 
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0.5 1.0 1.5 



b/b s 

FIG. 1: (Color online) Negativities between the first and the 
j th spin in an open spin s chain with first neighbor XY cou- 
pling, as a function of the transverse field b, with border cor- 
rections (see text), and three different values of s. We have 
set an anisotropy v y /v x = 1 — S/(2sn), with S = 2.5 and 
n = 8 spins. The factorizing field corresponds to the last 
parity transition, and is singled out as the field where all neg- 
ativities merge, approaching common non-zero distinct side 
limits. The lowest line (in red) depicts the end-to-end nega- 
tivity (jVi_„). 

before the last transition (roughly an W-state). b s plays 
in this small chain the role of a quantum critical field. 

Fig. [5] depicts results for a greater anisotropy S — 7.5. 
In this case just the last two parity transitions are visible 
in the negativity. The common side limits of N^j are 
smaller and all negativities exhibit a maximum to the 
right of the factorizing field. The behavior is then more 
similar to that of larger XY systems [HI]. Nonetheless, 
there is still a clear interval of full entanglement range 
around b s , with finite side- limits at b s when observed in 
detail (inset). 

The side limits at separability can actually be modi- 
fied in this system by changing the even-odd field ratio 
i] = b e /b , according to Eq. (|26[) . Results for a fixed ratio 
j] = 10 (with pertinent border corrections) are shown in 
Fig. [3J in which case separability is exactly attained at 
an odd field b os = b s /y/rj. We have again plotted just the 
negativities between the first and the j spin, which now 
approach two common side limits at each side, one for j 
even (N^ e ) and one for j odd (N^ ). While the former be- 
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FIG. 2: (Color online) Same details as fig. 1 for 5 — 7.5 and 
s = 3/2. The inset depicts the behavior in the vicinity of the 
separability field. 
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FIG. 3: (Color online) Same details as fig. 1 for 8 = 7.5, 
s = 3/2 and an alternating field with fixed even/odd ratio 
b e /b — 10. Red (blue) lines depict results for j even (odd). 



come quite small, the latter become clearly appreciable, 
the final effect for such large ratios being essentially that 
just odd sites become uniformly entangled in the vicin- 
ity of b os . Even-even negativities iVjj (not shown) are of 
course also very small at b os . Notice finally that N13 can 
become much larger than N12 in the region around b os , 
despite the absence of second neighbor couplings. 

In summary, we have first determined the conditions 
for the existence of separable parity breaking (and locally 
coherent) eigenstates in general XYZ arrays of arbitrary 
spins in a general transverse field, showing in particu- 
lar the possibility of exact separability in open as well 
as non-uniform chains through non-uniform transverse 
fields. We have also determined the entanglement prop- 
erties of the associated definite parity states, through 
the evaluation of the concurrence and negativity for any 
pair of spins or subsystems, for any spin values. These 
states, which approach both GHZ and VF-states in par- 
ticular limits, exhibit full entanglement range when non- 
orthogonal, and can be seen as effective two qubit entan- 
gled states for any bipartition. Moreover, the same holds 
for their uniform mixture as well as for the reduced den- 
sity of any subsystem. The finite entanglement limits at 
the SP become relevant in finite arrays close to the XX Z 
limit, where the separability field can be clearly identified 
with the last GS parity transition, as verified in the mi- 
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merical results presented, playing the role of a quantum 
critical field. The possibility of exact separability in an 
alternating field (b e = 7]b ) for arbitrary even-odd ratios 
rj, leading to controllable entanglement side-limits, has 
also been disclosed. The present results provide a deeper 



understanding of the behavior of pairwise entanglement 
in finite XYZ spin arrays subject to transverse fields. 
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CONICET (NC, JMM) of Argentina. 
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